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Abstract 

The 7 p —> j'qp reaction in the Sn (1535) resonance region is investigated as a method 
to access the 5'ii(1535) magnetic dipole moment. To study the feasibility, we per¬ 
form calculations of the yp —> 'jrjp process within an effective Lagrangian approach 
containing both the Sn resonant mechanism and a background of non-resonant con¬ 
tributions. Predictions are made for the forthcoming experiments. In particular, we 
focus on the sensitivity of cross sections and photon asymmetries to the S'ii(1535) 
magnetic dipole moment. 
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1 Introduction 


The study of the properties of the S'ii(1535) N* resonance, the lowest lying 
= \ nucleon resonance, provides valuable insight into the nature of QCD 
in the non-perturbative domain. In particular, the large mass splitting between 
the nucleon ground state iV(939) and its negative parity partner Al*(1535) is 
connected to the spontaneous breaking of the chiral symmetry. Indeed, if the 
two-flavor chiral symmetry were exact and preserved by the QCD vacuum, 
QCD would predict parity doublets degenerate in mass [1,2]. 

In recent years, new techniques have become available to perform numerical 
calculations of the properties of nucleon excited states using lattice QCD [3-5]. 
In particular, calculations of the mass spectrum of nucleon excited states were 
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performed [3], which implement chiral symmetry at non-zero lattice spacing, 
by use of domain wall fermions. The recent quenched lattice QCD calcula¬ 
tions [3-5], using different fermion actions, all reproduce quite well the ob¬ 
served large mass splitting between the nucleon and its negative parity part¬ 
ner iV*(1535). Furthermore, the next generation of lattice QCD calculations 
opens up the exciting prospect to study how nucleon and resonance proper¬ 
ties change when varying the quark mass parameter from the large quark mass 
regime, where chiral symmetry is explicitly broken, down to sufficiently small 
quark mass values, where one can use the chiral symmetry of QCD to perform 
the extrapolation to the chiral limit [6]. The large quark mass region may be 
the regime where a naive constituent quark picture could hold, and one may 
investigate how such a picture breaks down when approaching physical quark 
masses. 

A particularly useful observable to study in this respect is the baryon magnetic 
dipole moment (MDM), as it has often been questioned and discussed [7- 
9] why the naive non-relativistic constituent quark model yields a relatively 
good agreement with the experimentally known values for the MDMs of the 
ground state octet baryons. Recently, it has been studied [10,11] how the 
proton and neutron MDMs change when extrapolating from large up- and 
down-quark masses (about 6-10 times their physical mass values) where 
current lattice QCD calculations are performed, down to small quark masses, 
where the non-analytic behavior in the quark mass can be calculated from the 
chiral symmetry of QCD. In particular, it has been noticed in Ref. [11] that the 
apparent success of the non-relativistic constituent quark model in reproducing 
the observed ratio of proton to neutron MDMs at the physical quark masses 
is rather coincidental, supporting the hnding of Ref. [7]. Therefore, it may be 
very revealing to extend such a study to the MDMs of nucleon resonances, 
such as the A(1232) and the S'ii(1535). 

The A(1232) MDMs have been investigated on the lattice at rather large quark 
masses [12], and very recently the chiral extrapolation of the A(1232) MDMs, 
including the next-to-leading non-analytic variation with the quark mass, was 
also studied [13]. As there is a large gap in quark mass to bridge between the 
present quenched lattice QCD calculations and the chiral limit, it would be 
extremely helpful to know the resonance MDM for the physical quark mass 
values, through experiment. Unfortunately, the experimental information on 
the MDMs beyond the ground state baryon octet is very scarce. With the 
notable exception of the baryon, these higher nucleon resonances decay 
strongly, and thus have too short lifetimes to measure their MDMs in the 
conventional way through spin precession measurements. 

For the A’''+(1232), it was therefore proposed to measure its MDM through 
the 7r~^p 'jTi^p reaction [14], and two experiments were performed [15,16]. 
Using different theoretical analyses, the Particle Data Group (PDG) [17] 
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quotes as range of the MDM : /iA++ = 3.7 — 7.5 /xat, where is the nu¬ 
clear magneton. The large uncertainty in the extraction of /iA++ from the 
data is due to large non-resonant contributions to the 7r~^p 'yn~^p reaction 
because of bremsstrahlung from the charged pion (tt’*') and proton (p). Al¬ 
ternatively, it has been proposed [18] to determine the magnetic moment of 
the A’''(1232) through measurement of the 'jp reaction, thus avoiding 

bremsstrahlung contributions from the pions. First calculations for this reac¬ 
tion, including besides the A(1232) resonant mechanism [19,20] also a back¬ 
ground of non-resonant contributions, have recently been performed [21,22]. 
Due to the small cross sections of the 'jp —> reaction, the hrst experimen¬ 

tal data have been obtained only very recently by the A2/TAPS Collaboration 
at the MAMI accelerator [23]. The pioneering experiment of Ref. [23] sees a 
clear deviation from the known soft bremsstrahlung processes (for which the fi¬ 
nal photon is radiated from the external proton lines), and the measured cross 
sections are in qualitatively good agreement with the calculation of Ref. [21]. 
The deviation from the soft bremsstrahlung contains a sensitivity to the mag¬ 
netic dipole strength at the A(1232) resonance position. The experiment of 
Ref. [23] does not yet allow for a reliable extraction of this quantity, because of 
limited statistics. However, a dedicated experiment to measure the yp — 
reaction is planned at MAMI in the near future [24]. This experiment will use 
the dTT Crystal Ball detector and thus increase the count rates by two orders of 
magnitude. In addition, this forthcoming experiment [24] will be performed 
with a polarized photon beam, because the sensitivity of the yp —yvr^p reac¬ 
tion to the A’''(1232) MDM can be increased substantially when using linearly 
polarized photons, as has been shown in Ref. [21]. 

In the present paper, we explore the extension of this technique to access the 
MDM of the S'ii(1535) resonance. Since the S'ii(1535) resonance is known 
to dominate the yp —pp reaction in the threshold region, the yp —ypp 
reaction, with a low energy outgoing photon (up to about 100 MeV), seems to 
be a promising tool to isolate the MDM of the An (1535) resonance. Such an 
experiment can be envisaged using the Crystal Barrel@ELSA and in the near 
future with the Crystal Ball@MAMI, once the higher beam energy of about 
1.5 GeV becomes available at MAMI-C. 

In Section 2, we provide some model estimates for the An(1535) MDM. In 
particular, we present the calculation of the An (1535) MDM in the constituent 
quark model and contrast it with the value taken in the picture of the An (1535) 
as a {KTi) quasi-bound (molecular) state, as has been proposed in Ref. [25]. 

We then develop in Section 3 an effective Lagrangian model for the yp —ypp 
reaction in the An (1535) region, which is obtained by coupling a photon in 
a gauge invariant way to an analogous model for the yp —pp reaction. The 
model contains both the An (1535) resonance contribution and a background 
of non-resonant processes. 
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In Section 4, we discuss the model-independent link between the 7 p —rjp 
process and the 7 p —'yrjp process in the limit where the outgoing photon 
energy approaches zero. We present a low energy theorem, which provides a 
strong constraint for both the experiment and the theoretical description of 
the 7 p — 7 ? 7 p reaction. 

In Section 5, we then use the developed effective Lagrangian model as a tool 
to study the feasibility of an experiment to measure the jp 'ypp reaction. 
In particular, we investigate the sensitivity of the 7 p —'jpp cross sections and 
photon asymmetries to the S'ii(1535) MDM. The hnal section summarizes our 
hndings. 


2 Model calculations for the S'ii(1535) magnetic moment 


In this section, we calculate the results for the MDM of the S'ii(1535) in 
the framework of the constituent quark model, and within the picture where 
the S'ii(1535) is dynamically generated as a quasi-bound state, as hrst 

proposed in Ref. [25]. 


2.1 Magnetic moments in the constituent quark model 


In the nonrelativistic SU(6) constituent quark model, the lowest-lying neg¬ 
ative-parity nucleon resonances are \N‘^Pi/ 2 ) and \N'^Pi/ 2 ), where the usual 
spectroscopic notations ^Pi /2 and ^Pi /2 are used to indicate their total quark 
spin S = 1/2, 3/2 {2S -|- 1 = 2, 4), orbital angular momentum L = 1 {P- 
wave), and total angular momentum J = 1/2. In contrast, the ground-state 
baryon octet (e.g., N, A, S) and decuplet (e.g.. A) have the same S'-wave 
spatial wavefunctions with L = 0, and thus J = S. The wavefunctions of the 
l^^-Pi/ 2 ) and lA^Pi/ 2 ) states are given explicitly as 


\n^Pi/2) = ^ E ( 11 1H ) 


mirris 


X { 


72 (Am. 0 '’ + Xm. 0 ^ 


+ tlj 


Imi 


71 0 '’- 0 ^) 


( 1 ) 


\N^Pi/2) = E ( 1 I l| I) 01m, Am. + 0L, Am. 0 

V ^ mirris 


( 2 ) 


where y, and 0 denote the spatial, spin, and flavor wavefunctions. The 
superscripts s or p (A) of these wavefunctions indicate that they are totally 
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symmetric among three quarks, or odd (even) under the exchange of the first 
two quarks. Further details can be found in Ref. [26]. 

However, the observed lowest-lying negative-parity nucleon resonances are the 
S'ii(1535) and 5'ii(1650), obtained as conhguration mixtures of the liV^Pi/ 2 ) 
and liV^Pi/ 2 ) SU(6) states. 


1 All(1535)) = 1 N‘^Pi/2 ) cost? - I N^Pi/2 ) siuT? , 

1 All(1650) ) = 1 N‘^Pi/2 ) siuT? + I N^Pi/2 ) cost? , 

( 3 ) 

where t? denotes the mixing angle. 


In a constituent quark model, the magnetic dipole moments of qqq baryons 
consist of contributions from both quark spin and orbital angular momentum, 
i.e. /X = with 

^ ^ nrii 

( 4 ) 

i i 

( 5 ) 


where the index i sums over three quarks. 

In Appendix A, we present the calculation of the spin and orbital contributions 


to the magnetic dipole moments for the states iV^Pi/2) and iV^Pi/2) with 
Jz = - 1 - 1 / 2 . The results are given by : 

/i(?V'Pl/2 + ) = ( iV2Pi/2+ 1 liz 1 iV2Pi/2+ ) 

g (2/Tt( -|- fJ^d) 3 j 

/i(?V'Pi/2 °) = ( N^Pi/2 ° 1 /i. 1 N^Pi/2 ° ) 

( 6 ) 

9 (2/T(T “|- fJju) 0 ) 

fi{N^Pl/ 2 ^) = ( iV"Pl/2+ 1 fiz 1 N^Pi/ 2 ^ ) 

( 7 ) 

9 T h'rf) 3 h'V ) 

/i(iV"Pi/2 °) = ( N*Pi /2 ° 1 /i. 1 N*Pi /2 ° ) 

(8) 

9 T fd'u) 3 fd'N ) 

( 9 ) 

where /r.^ = /rf -|- /if', and the superscripts -|- 

and 0 denote the charge of 


the resonance state. Here we assume the constituent quarks as Dirac point 
particles and approximated = rud = \ miq. Using the relation jj, = Q/2m, 
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we obtain Hu = Qul^^'rnu = 2 /xat and Hd = Qd/‘^md = —fiN, which are used to 
obtain the values in Eqs. ( 6 )-(9). 

In addition, there are the cross terms for /if which are obtained from mixing 
the |iV^Pi/ 2 ) and |iV^Pi/ 2 ) states, 


( 7V^Pi/ 2+ I /xf I N^Pi/2^ ) = { N^Pi/2^ I /xf I N^Pi/2^ ) 

9 (/^w f^d) fd-N 7 (10) 

( N^P,/2 ° I /if I N^Pi/2 ° ) = ( N^Pi/2 ° I /if I N^Pi/2 ° ) 

9 {.f^d f^u) 3 fd'N • (11) 

On the other hand, there are no cross terms for /if because \N‘^Pi/ 2 ) and 
\N'^Pi/ 2 ) have orthogonal quark spin states which are not affected by /if. 

Now the magnetic moments of the S'ii(1535) and S'ii(1650) resonances can 
be expressed in terms of the magnetic moments of the |iV^Pi/ 2 ) and |iV^Pi/ 2 ) 
states and the cross terms. 


/i(S'ii(1535)) = /i(iV^Pi/2) cos^ d + /i(iV^Pi/2) sin^ d (12) 

- 2 ( iV^Pi /2 I /if I iV'^Pi /2 ) sin cos d , 

/i(^ii(1650)) = /i(iV2Pi/2) sin^ d + /i(iV^Pi/2) cos^ d (13) 

+ 2 ( iV^Pi /2 I /if I fV"‘Pi /2 ) sin d cos d . 


The value of the mixing angle d depends on the quark interaction. Assuming 
a hyperhne interaction between the quarks, Isgur and Karl [26] predicted a 
mixing angle -d = tan“^ (v^— l)/2 — —31.7°, which is close to the empirical 
mixing angle d ~ —32° found in Ref. [27]. Using the value d = —31.7° in 
Eqs. (12) and (13), we obtain 


/^ 5 +pi 535 ) — 1-89 /iAT and /i 50 ^( 1535 ) — —1.28 /xat , 

/^5+(i650) = 0.11 /iAT and /isopieso) = 0.95 , 

which agrees with the result by the Genova group [28]. Very similar results 
can also be obtained when the interactions of the S'!! resonances coupling to 
the 71N and r]N channels are explicitly considered [29]. 
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2.2 Magnetic moment o/S'ii(1535) as a quasi-bound state 


It has been proposed in Ref. [25] that the S'ii(1535) can be considered as a 
(RE) quasi-bound state of mesons and baryons. Using an SU(3) effective chiral 
Lagrangian and iterating the derived coupled-channel potentials to inhnite 
orders, the S'ii(1535) resonance can be generated dynamically. This hnding 
has also been verihed within the chiral unitary approach [31,32]. Especially, 
in Ref. [32] it is found that the coupling of the S'ii(1535) to the rjN channel 
is nearly as large as the one to the KTi channel. 

In a simplihed picture, we regard the S'ii(1535) as a (R"S) quasi-bound state. 
Then the S'ii(1535) charge states can be written as : 


S+(1535) > = (A'+ E») + (A-» E+), 

S;i( 1535)> = y|(/f+E-) - (14) 

From Eq. (14), it then follows that the magnetic moment of both 5'ii(1535) 
charge states can be expressed in terms of the hyperon magnetic moments. 


1 2 

^■5+(1535) “ 2 A'-SO + 2 h-Sl- , 

2 1 , , 

hSh(1535) — 2 ^^” 3^^°‘ 

Using the experimental values for the S'*' and magnetic moments, and 
approximating the unknown magnetic moment by the average of the S’*' 
and magnetic moments, one obtains for the 5'ii(1535) magnetic moments 
in the (R"E) molecule picture the values : 


^■5^(1535) — 1-86 /iAT , 

/^Sb(1535) = -0.56 /iiv. (16) 


We like to stress that the calculation given here for the MDM of the 5'ii(1535) 
as a (R"E) quasi-bound state is only a very simple exercise of the MDM for 
a dynamically generated resonance, and a full calculation within a coupled- 
channel framework remains to be done. In this respect, the magnetic moment 
of the A(1405), which is in the same lowest energy = \ baryon octet as 
the All (1535), has recently been evaluated in Ref. [33] using unitarized chiral 
perturbation theory. In this approach the A(1405) is dynamically generated. 
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It would be very worthwhile to also perform such a calculation for the MDM 
of the S'ii(1535) resonance. 


3 Effective Lagrangian model 


In this section we develop an effective Lagrangian model for the yp —'yrjp 
reaction in the 5'ii(1535) resonance region, which will subsequently be used 
as a tool to investigate the sensitivity to the S'ii(1535) MDM. 

In the 7 p —ypp process, a photon {k, A) hits a proton target {pn, sn)', in 
the hnal state, a photon {k', A'), an p meson (g), and a proton (p'^, s'j^) are 
observed. Here k, k', p^, p'j^, and q are the four-momentum for the respective 
particles, A and A' denote the photon helicity, and sn and s'^ are the proton 
spin projection. 

For convenience, our results for the experimental observables will be expressed 
in the center-of-mass (c.m.) frame of the initial yp system. The kinematics of 
the yp —ypp reaction can be described by hve variables. First, we choose 
the energies and of the initial and outgoing photon, respectively. 
The other three variables are the polar and azimuthal angles 6^'^, of the 
hnal photon, and 6^'^ the polar angle of the eta meson. They are dehned by 
choosing the plane which contains the initial particles and the hnal p to be 
the plane with k pointing in the 2 ;-direction, and thus = 0. 

The unpolarized hve-fold diherential cross section for the yp —ypp reaction, 
which is diherential with respect to the outgoing photon energy and angles as 
well as the eta meson angles in the c.m. system, is given by 

da _ 1 1 |q|2 

dE;"” dD™ “ (27r)5 32v^ | q | (E; + Ug) + E'^'^Ug cos 

X I 7 E E E EI C(k', A') Al-'- M , (17) 

where lOg and q are the energy and momentum of the eta meson, E'^ is the 
hnal proton energy, the c.m. angle between the outgoing photon and the 
eta, and e^{k^ A) and e*y{k'^ A') are the polarization vectors of the incoming and 
outgoing photons, respectively. Furthermore, is a tensor for the yp —a ypp 
process, which will be calculated next within an ehective Lagrangian model. 

To develop a model for the yp —a ypp reaction, we start from an ehective 
Lagrangian model for the p photoproduction reaction yp —a pp, whose param¬ 
eters are determined by a ht to the experimental data. The hxed parameters 
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Fig. 1. Diagrams for the 7 p —> r/p reaction in the 5ii(1535) region : 5ii resonance 
excitation (a), vector meson exchange (b), and Born diagrams (cl,c2). 

are subsequently used to calculate the yp —ypp process. The dominant con¬ 
tributions for Tj photoproduction in the S'ii(1535) region are shown in Fig. 1. 
These include, in addition to the resonance contribution, the nucleon Born 
terms and the t-channel vector meson exchanges as the nonresonant back¬ 
ground. Besides the dominant 511(1535) nucleon resonance, we also include 
the 5ii(1650) since its contribution in the 5ii(1535) region is not negligible. 
Other nucleon resonances are neglected at this stage since their contributions 
to the r] photoproduction cross sections in the 5ii(1535) region are small. 


To extend our description to the yp —>• ypp reaction, we couple a photon 
in a gauge invariant way to all charged particles in Fig. 1. The resulting dia¬ 
grams are shown in Fig. 2, where the 7-5ii-5ii vertex in diagram (a2) contains 
contribution from the magnetic moments of the 5ii resonances. Here we as¬ 
sume that the anomalous MDM of the 5^(1650) vanishes, because it has only 
a very small effect in the 5ii(1535) region and its predicted value from the 
quark model is small (see Sec. 2.1). Therefore, in comparison with the yp —>■ pp 
process, the only new parameter entering in the description of the yp —>■ ypp 
process is the 5^(1535) anomalous magnetic moment i^s+{ 1535 )- follow¬ 

ing, we discuss the 5ii resonance as well as the background contributions to 
the yp —> ypp reaction. 


3.1 Sii{1535) resonance propagator 


As the 5ii(1535) resonance is a spin-1/2 particle, its Feynman propagator is 
given by : 

GH(p) = j-d—, ( 18 ) 

p- Mr 

where p is the four-momentum and Mr the mass of the 5ii(1535). To take 
account of the hnite width of the 5ii resonance, we follow the procedure of 
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Fig. 2. Diagrams considered in the calculation of the jp —> j'qp reaction in the 
5ii(1535) region, obtained by gauge invariant coupling of a photon to the diagrams 
of Fig. 1 : Sii resonance diagrams (al-a3), vector-meson exchange diagrams (bl-b2), 
and Born diagrams (cl-c6). 

Refs. [34,35] by using a complex pole description for the resonance excitation. 
This amounts to the replacement 

Mr ->■ Mr — - Fr (19) 

in the propagator of Eq. (18). This ‘complex mass scheme’ guarantees electro¬ 
magnetic gauge invariance. In contrast, the use of a Breit-Wigner propagator 
with an energy-dependent width will violate gauge invariance when applied to 
the Sii contribution for the yp —> yyp reaction. By fitting the yp —> rjp data 
in the 511(1535) region, we can determine the complex pole of the 5ii(1535). 
The £t results are shown in Fig. 3. Limiting our £t to the data near the 
5ii(1535) peak in the range = 730 — 880 MeV, we obtain [Mr, Tr) = 
(1490, 154) MeV, which will be used in our calculation for the yp —> yyp 
process. In comparison, the value given by the PDG [17] is (1505, 170) MeV. 
However, using the PDG value the agreement with the photoproduction data 
is less satisfactory. 


3.2 5ii(1535) resonance amplitudes 


The interaction Lagrangians relevant for the resonance contribution {R = Rn) 
in the yp ^ yyp process are 
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Fig. 3. Total cross section for the 'jp rjp. The solid line is the full result, and 
the dashed line indicates the contribution from the S'ii(1535). The data are from 
TAPS [36] and GRAAL [37], 


= -ei)R[ ^/JR , 

^yNR = - 2(m/+Mk) ^NR 75 + h.c., 

^^■nNR = — i g-^NR 'ipR 0r? + h.C. , 


( 20 ) 


where g^^NR is the pseudoscalar coupling for the gNR vertex, cr and kr the 
resonance charge and the anomalous magnetic dipole moment, krir the N 
R transition magnetic coupling, and — d'^A^^. 


Using these effective Lagrangians, we can derive the amplitudes for the S'!! 
resonance contributions in Fig. 1(a) and Fig. 2(al-a3). The amplitude for the 
S'ii(1535) resonance contribution to the 7 p —gp process, corresponding to 
Fig. 1(a), is given by 


e^{k,\)M’p{gp gp) 

= i s)v) Gr{pr) KpRg^ N{pn, sn) , (21) 

where Gr is the S'!! propagator given by Eqs. (18) and (19). 

Starting from Eq. (21) to the gp gp process, we now construct the corre¬ 
sponding Fii contribution to the gp ggp process. This is obtained by cou¬ 
pling a photon in a gauge invariant way to all charged particles in Fig. 1(a). 
First we obtain the diagrams with a photon attached to an external proton 
[Figs. 2(al) and 2(a3)], described by the usual Dirac and Pauli currents with 
Kp = 1.79 for the proton anomalous MDM. The contributions of these dia¬ 
grams to the tensor of Eq. (17) are given by 
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( 22 ) 


-^al(7P^ 7^p) 


= i ka N{pj^, s)v) Gnipj^) KpR 75 a'"" 


Mn+Mji 

X 5'jv(pAr - k') 

MZilP^lVV) 


* Mjv+Mfl N{Pn^ ^n) 




2M^ 


N{pn,sn) , 


7 S 2Mjv 


X Sn{p'n + k') Gr{pr) KpR 75 a'"" iV(pAr, Sjv) , 


(23) 


where p'p, and pr correspond to the four-momenta of the S'!! resonance in 
Fig. 2(al) and 2(a3). In the soft-photon limit {k' 0), the coupling of the 

photon to the external lines gives the only contribution. However, at hnite 
energy for the emitted photon, gauge invariance also requires the diagram 
with a photon attached to the intermediate S'!! resonance [Fig. 2(a2)], which 
corresponds to the tensor. 


MlYw IPP) 

= * wSm-a "'iv) Gr{pY T^rr Gr{pr) 

X KpR 75 0-^“ N{pn: Sn) , ( 24 ) 

with Gr{pr) and Gr{p'Y) the S'!! propagators before and after emitting the 
photon in Fig. 2(a2), and the 'jRR vertex 

Prs = 7" - i'tKtr'"’Jj, (25) 

containing the photon coupling to the charge of a spin 1/2 held of a Dirac 
particle [34] as well as the Pauli current proportional to kr. 

The vertex for a point spin-1/2 particle follows from the requirement of gauge 
invariance. Indeed, only the sum of the three diagrams (al, a2, a3) in Fig. 2 is 
gauge-invariant, which is expressed through the electromagnetic Ward identity 
relating the ■yRR vertex with the propagator, 

{Pr - Pr)u G{Pr) G{pr) = G{pr) - G{pY> . (26) 

To include the hnite width of the 5'ii, the ‘complex mass scheme’ was advo¬ 
cated in Ref. [34], i.e., the replacement of Mr by the complex mass of Eq. (19). 
In doing so, it is easily checked that the Ward identity of Eq. (26) still holds. 
On the other hand, when using Breit-Wigner propagators, and replacing 
in the denominator of Eq. (18) by M|, — i MrTr(W), with Fj:j(IF) an energy 
dependent width, one immediately sees that the propagators before and af¬ 
ter the emission of the photon have diherent widths (except for the soft-photon 
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limit k' 0, where the contribution of diagram Fig. 2(a2) vanishes relative to 
those of Fig. 2(al) and (a3)). Therefore, when using energy-dependent widths 
at hnite energy of the emitted photon, it is not possible to maintain the Ward 
identity of Eq. (26) with the vertex of Eq. (25). 


3.3 Background contributions 


The effective Lagrangians relevant to the background contribution in the yp 
'jrjp process are 


k^'rNN = -ei^N 'ipN , 

k^r)NN = —i drjNN '4’N 75 '4’N 4>v > 

^ W 4>, {d^V-) , 

C-VNN =gVNN'ipN ’ 


(27) 


where V denotes p and ui vector meson helds, and the pseudoscalar coupling 
QrfNN is used in our calculation. With these effective Lagrangians, it is straight¬ 
forward to derive the amplitudes corresponding to the diagrams of Fig. 1(b) 
and (cl,c2) as well as Fig. 2(bl,b2) and (cl-c6). 


For the pNN coupling in Eq. (27), we determine its value by fitting the photo¬ 
production data and obtain = 0.05. For the vector meson couplings, 

we use the values from Ref. [38]. There are uncertainties associated with the 
grjNN and vector meson couplings, and it could affect the correct extraction of 
the S'ii(1535) MDM. However, the background contributions, which contain 
the Born terms and vector meson exchanges, do not exceed 1% of the yp —pp 
and yp —ypp cross sections in the S'ii(1535) energy region considered here. 
Therefore, we believe that these uncertainties will not signihcantly affect the 
extraction of the S'ii(1535) MDM. 


4 Soft photon limit 


In the soft photon limit {k' —> 0), gauge invariance provides a model-inde¬ 
pendent relation between the cross sections for the yp —ypp and yp —pp 
processes. In this section, we formulate this low energy theorem. The expres¬ 
sion we derive applies to the production of any neutral meson M, and we 
therefore refer more generally to the yp —yMp reaction in this section. In 
particular, it also applies to the yp —y7r°p reaction as studied in Ref. [21]. 
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Fig. 4. Angular weight-function W(v) of Eq. (29) in the soft-photon formula for the 
7P IVP (solid curve) and 7 p ^ 77r°p (dashed curve) processes, as function of the 
meson angle 0^. 

All variables in this section are expressed in the c.m. frame of the initial 'jp sys¬ 
tem. In the soft photon limit, the momentum and energy of the outgoing pho¬ 
ton vanishes, and the jp 'yMp reaction is dominated by the bremsstrahlung 
process from the initial and hnal protons. In this limit, when integrating the 
hve-fold differential cross sections for the 7 p —'yMp process over the outgoing 
photon angles, we obtain (the details are given in Appendix B) 


Ml _. J-ILwm (—) 


as /c' —0 . 


(28) 


Here {da/dVluJ^ is the differential cross section for the 7 p —Mp process, 
and we dehne an angular weight-function 


W{v) = -1 + 



In 


n -I- 1 


(29) 


where v = Jl + with t = {k — qY and Mjv the nucleon mass. 


If we integrate the cross section of Eq. (28) also over the meson angles, the 
result is 


= f dPM 

dE'^'^ J ^ dE'^'^dklf^ ’ 

—^ — Mm as /c' ^ 0 , 

^/cm ’ 

7 

with a “weight-averaged” total cross section, 

2 

a« = ^/ da'sw(v) 



(30) 


(31) 
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Fig. 5. Outgoing photon energy dependence of the ratio R, defined in Eq. (32), for 
the 7 p —> yr/p reaction at = 900 MeV. 

We show the angular weight-function W{v) in Fig. 4 for both the 7 p —>■ 7 ? 7 p 
and the 7 p ^ 'jir^p processes. If the meson is emitted in the forward direc¬ 
tion with respect to the initial photon direction, there appears a destructive 
interference between the bremsstrahlung from the initial and hnal nucleons, 
whereas for a meson emitted in the backward direction both bremsstrahlung 
sub-processes add constructively. 

The low energy theorem of Eq. (30) provides a check for both theoretical 
model calculations and experimental measurements, because the ratio 

R = ■ {daa (32) 

with da/dE'df'^ for 7 p —> yMp and dM for 7 p —>■ Mp, approaches unity when 
k' 0. 

In Fig. 5, we show the ratio R constructed from our theoretical calculations of 
the 7 p —>• 'jpp and 'jp —> pp reactions. We indeed observe that R approaches 
unity when k' 0. For k' ~ 20 MeV, this ratio is still within 5% of the soft 
photon limit. Upon increasing the outgoing photon energy further, one starts 
to see clear deviations (of the order of 30%) from the soft photon limit. These 
deviations allow one to extract new resonance information from the jp —> ppp 
reaction, which will be studied in the next section. 


5 Results for 'yp —> ^pp observables 


In this section, we use our effective Lagrangian model as a tool to study the 
feasibility of an experiment to measure the pp 'ypp reaction. In particular, 
we present our results for the 'yp —>■ 'ypp cross sections and photon asymmetries 
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Fig. 6. Outgoing photon energy dependence of the cross section da/dE'^™ (inte¬ 
grated over photon and eta angles) for the yp ^ yr/p reaction at = 900 MeV. 

for different values of in the range of the estimates made in Section 2, 

i.e., between —2 to +2 (in units of ptv), to illustrate the sensitivity with respect 
to the S'ii(1535) MDM. 

In Fig. 6, we show the outgoing photon energy dependence of the cross sec¬ 
tion da/dE'^'^ for the yp —ypp reaction integrated over the photon and 
eta angles for fixed = 900 MeV. The cross section at low energies E'^'^ 
exhibits the characteristic bremsstrahlung behavior, but displays clear devia¬ 
tions from such a behavior when going to higher outgoing photon energies (see 
also Fig. 5). To give an idea about the feasibility to measure the yp —ypp 
cross section, it is worthwhile to make the comparison with the yp —yvr^p 
reaction in the A(1232) region, which has already been measured [23]. At zero 
outgoing photon energy, the cross section for radiative eta production can be 
estimated from the yp —pp cross section. On the other hand, the maximum 
value of the yp —pp cross section at the S'ii(1535) position is about 5 % of 
the yp —7r°p cross section at the A(1232) position. However, at an outgoing 
photon energy of about 100 MeV, which is the region where it is sensitive 
to the MDM of the resonance, our estimate for the yp —ypp reaction at 
E^^^ = 900 MeV yields count rates that are about 10-15 % of the measured 
value of Ref. [23] in case of the yp —yTT^p reaction at E]^^ = 450 MeV, which 
turns out to be much more favorable than the 5 % found in the soft photon 
limit. 

One can increase the sensitivity of the yp — ypp reaction to the An (1535) 
MDM by measuring differential cross sections in specihc kinematical regions. 
In Fig. 7, we show the cross section, in the c.m. frame of the initial yp system, 
for different values of the p angle upon integration over all photon directions 
at E^^'^ = 900 MeV. At the lower energies, one sees a strong forward-backward 
asymmetry, which was discussed before in Fig. 4, yielding substantially larger 
cross sections for backward p c.m. angles. However, at E'd^'^ ^ 100 MeV, the 
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Fig. 7. Outgoing photon energy dependence of the . d?a /for the 

IP IVP reaction at = 900 MeV. 

cross sections become comparable for all angles and exhibit a 20 % change 
when Ksjpisss) is varied between —2 to +2. 

One can fnrther increase the sensitivity by measnring cross sections at par- 
ticnlar angnlar regions of the emitted rj and 7 . As an example, we show in 
Fig. 8 the hve-fold c.m. differential cross section d^a /and the 
photon asymmetry S at = 900 MeV for the polar angles = 90° and 
^cm _ three valnes of the azimnthal angle = 0°, 90°, and 180° 

(note that the rj is emitted in the plane, i.e., = 0°). We notice that the 

photon asymmetries reach valnes of abont 10% aronnd E'^'^ = 100 MeV, and 
show a strong dependence on the ont-of-plane angle of the emitted photon. 
One also hnds that the sensitivity of the photon asymmetry to the valne of 
^Sii(i 535 ) is qnite sizeable. 


6 Summary and conclusions 


In this paper we investigated the jp 'ypp reaction in the An (1535) resonance 
region as a tool to access the An (1535) magnetic dipole moment (MDM). We 
performed calcnlations of the yp —ypp process nsing an effective Lagrangian 
approach, which contains the An (1535) and An (1650) resonant mechanism as 
well as a backgronnd of non-resonant contribntions. As a starting point, the 
effective Lagrangian formalism is hrst applied to the yp —pp reaction. When 
extending the calcnlation to the yp —ypp process, a photon is conpled to 
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E’™ [MeV] [MeV] E’^"" [MeV] 

Fig. 8. Outgoing photon energy dependence of and pho¬ 

ton asymmetry S of the yp ^ reaction for fixed polar angles 0™ = 90° and 
^cm _ varying azimuthal angle (p'^ = 0°, 90°, 180° (from left to right) at 

E}^^ = 900 MeV. 

all charged particles in the jp pp process, and the coupling is constrained 
by gauge invariance. In particular, to take account of the finite width of the 
S'!! resonance, a complex mass scheme rather than a Breit-Wigner propagator 
with energy-dependent widths, has been used to maintain the Ward identity. 

We presented our results for the pp ppp process using both unpolarized 
and polarized photon beams. In particular, we focused on the sensitivity of 
cross sections and photon asymmetries to the S'ii(1535) anomalous MDM. We 
derived a low-energy theorem relating the cross sections for the pp pp and 
IP IVP processes in the soft-photon limit, which provides a useful check for 
both theoretical calculations and experimental measurements. 

We hrst showed the unpolarized cross section integrated over the full pho¬ 
ton and p angles. At low energies this cross section is dominated by the 
bremsstrahlung contribution, which is consistent with our check in the soft- 
photon limit. The sensitivity of the unpolarized integrated cross section to 
^Sii(i 535 ) is found to be small, but can be increased by measuring differential 
cross sections at particular angular ranges. For the differential cross section 
d?a /integrated over the photon angles (but not the p angles), and 
the hve-fold differential cross section dP’a/dE'^’^dVE^dVtp^, we found that the 
change when varying K,Siiii535) between —2 to +2 is at the 20%-25% level. One 
can further increase the sensitivity by measuring cross sections using a linearly 
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polarized photon beam. Indeed, our results clearly show that the sensitivity 
of the photon asymmetry to ^ 511 ( 1535 ) is quite sizeable. 

One inspiring motivation to study the MDMs of the S'ii(1535) resonance is 
to clarify the nature of this resonance—whether it is a traditional qqq state, 
or dynamically generated from meson-baryon interactions. In the framework 
of the SU( 6 ) quark model, the 5'ii(1535) and S'ii(1650) resonances are con- 
hguration mixtures of two SU( 6 ) states with excited orbital wavefunctions. 
Calculating both the quark spin and orbital angular momentum contribution 
for the MDM, as well as the cross terms due to the conhguration mixing, 
we obtained the MDM values for the 5'ii(1535) resonance. In comparison, we 
also presented a simple MDM calculation, where the 5'ii(1535) is dynamically 
generated and described as a molecule. Of course, further predictions 

of S'ii(1535) MDMs from other baryon resonance calculations would be most 
welcome, since through this reaction we can investigate the MDMs of orbitally 
excited nucleon resonances for the first time by an experiment. Although chal¬ 
lenging, such investigations seem now feasible at ELSA and MAMI. These 
planned experiments will provide an interesting test for every baryon struc¬ 
ture calculation, and give insight into the nature of the parity partner of the 
nucleon. 
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A Quark model calculation 


To perform a calculation in a constituent quark model with harmonic wave 
functions, it is convenient to work with Jacobi coordinates, which are dehned 
in terms of the three quark positions r*. 
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R=^{ri + r2 + Vs), 

P=^ir^-r2), 

A=^(ri + r2-2r3). 


(A.l) 


Their inversion relations are 


rs = R - ^X. 

We usually calculate in the center-of-mass system of the three quarks in a 
baryon, i.e., -R = 0. 

In a constituent quark model, the MDM of a baryon consists of contributions 
from both quark spin and orbital angular momentum, 



(A.3) 

Y i^i 

(A.4) 

y^p\ = y h, 

(A.5) 

Qi 

> - 'f'i X Pi, 

Qi + Q2 1 /Qi + Q 2 1 Qs ^ 

4m ^ \ 12m 3m'y 

+ /- (p X Px + Xx PpV 

4v3m ^ ^ 


where Qi denotes the charge of the ith quark, the index i sums over three 
quarks, and Ip = p x Pp, lx = X x p\. 

The calculation for the quark spin part of baryon MDMs is straightforward. 
But the calculation for the orbital angular momentum part is much less fa¬ 
miliar, since it vanishes for ground state octet and decuplet baryons, which 
have no orbital excitation. To calculate the matrix elements of with wave- 
functions of Eq. (1), we need to operate on spatial wavefunctions with 
2 :-components of l\, Ip, p x px, and A x Pp. We list in the following the 
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properties where jj,^ acts on spatial wavefunctions t/’im, 


{p X Px)zi^lm = 


(A X Pp),'tl) 1 ^ = mij 


Im ") 


^pz '4^1m 0 7 

hz Pirn = 

(p X Px)zPlm = 0, 

(A X Pp)zPL = 0- 


(A.6) 


In the case of a symmetric qqq state | {qqq)s), the calculation of p^ can be 
simplihed as 


Qg 

{qqq)s \ p \ {qqq)s) = ( {qqq)s \ — | {qqq)s) ■ 

ms 


(A.7) 


In the Tables A.1-A.4, we list the matrix elements for quark spin and orbital 
angular momentum contributions of the MDMs for the states | N‘^Pi/ 2 ) and 

I N‘^Pi /2 ) with Jz = +1/2. 

Table A.l 

Matrix elements for ( \ pf | +^^‘^"''^+ 1 / 2 ''") 


hf 

\N^Pi/2^) \N^Pi/2^) 

(iV'7^1/2 +1 
(iV"7’i/2+l 

~Q (2pti + Pd) 9 {Pu ~ Pd) 

9 {Pu ~ pd) 9 (2pn + Pd) 


Table A.2 

Matrix elements for ( I pf I N'^^~^^Pi/ 2 ^ ) 


kz 

\N^Pi/ 2 ^) \N^Pi/ 2 ^) 

{N^Pi/ 2 ^\ 

{N^Pi/ 2 ^\ 

~9 + fJ^u) 9 (/^d ~ f^u) 

9 (/^d l^u) 9 (2/dd H“ /^ti) 


Table A.3 

Matrix elements for ( N‘^^~^^Pi/ 2 ~^ \ Pz I ) 


kz 

\N^Pi/2^) \N^Pi/2^) 

{N^Pi/2^\ 

{N*Pi/2^\ 

K2+ + S) 0 

0 —k(zr + ^Si) 

J-O \ ''i'll ''^d ' 


B Derivation of the low energy theorem relating the 'yp —^ 'yMp 
and 'yp —Mp processes 


In the limit k' ^ 0, the 'yp 'yMp reaction is exactly described by the 
bremsstrahlung process from the initial and hnal protons. This yields the hve- 
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Table A.4 

Matrix elements for ( I I ) 


dz 

\N^Pi/2^) \N^Pi/2^) 

{N^Pi/2Y 

5 ( 2 ^ + !;) 0 

n 1 ( Qd 9 Qn \ 

^ 18 V md ^ m„ / 


fold differential c.m. cross section 


da 


/ da Y” 

^ p/cm 

■ £{k', A^) 

Pn ■ e{k', Xp 

\dQM J 

16M 

p'n ■ k' 

Pn ■ k' 


(B^l) 


as k' —0. In Eq. (B.l), {da/dkluY™ is the c.m. cross section for the qp —> Mp 
process, A..y = ±1 is the soft photon polarization, and e its polarization vector. 
We calculate the RHS of Eq. (B.l), by performing the sum over the photon 
polarizations, and integrating over the photon angles. This gives the result : 


da f da 


^ jpfcm j 

W ^ ' ' 


where we introduced the photon angular integral I as : 


I = 



‘^Pn-p'n 
{Pn ■ k') {pY ■ k') 


M 2 

{Pn ■ k'f 


M 2 

{p'n ■ k'f 


(B.2) 


(B.3) 


In Eq. (B.3), the second and third terms arise from the contribution of brems- 
strahlung from the same proton, whereas the hrst term stems from the inter¬ 
ference between the bremsstrahlung from the initial and hnal protons. 


We next work out the photon angular integral of Eq. (B.3). It turns out to be 
useful to introduce the initial and hnal nucleon velocities I3n = Pn/En and 
/3)y = p'j^/EY, and a Feynman parametrization in the hrst term of Eq. (B.3), 
which brings the two propagator factors to the same denominator. This leads 
to the expression : 


I = 


271 


(^/cm)2 


(1 — (3^ ■ /3^v) 


r+i 


1-1 


dy 


dx _ ^ _ 

1-1 {l-PyXY 


(B.4) 


(1 ~ Xk) Y 




dx ■ 


{1-YnxY ^ {l-YYxY}^ 

where Pn, P'n) Py are the magnitudes of Pn, P'n, and Py, which are 
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related by 


= /^AT -(1 + y) + f3'N-{l-y). (B.5) 

The x-integrals in Eq. (B.4) can be worked out using 

which yields for Eq. (B.4) the result : 

^ ^ ■ {-4 + 2 (1 - / 3 . . / 3 W ^ ( B . 7 ) 

By dehning the variable v = with t = {k —puY — (Pn ~PnY, we 

express 

(1 - / 3 iv ■ / 3 W = (1 - • (1 - . ( B . 8 ) 


and after a little algebra we can work out the remaining integral using : 


r+l 


dy 


-1 "ii-Pi) ii - 


2v 


2 In 


n + 1 


(B.9) 

Combining Eqs. (B.7-B.9), we obtain the result for the photon angular integral 


I = 


271 


^E'cray 


■4-i-l + 


v^ + l 
2v 


In 


V + 1 


(B.IO) 


Using Eq. (B.IO), the soft photon limit for the cross section of Eq. (B.2) then 
leads to the result of Eq. (28). 
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